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Abstract 

We analyze the quasinormal modes of L>-dimensional Schwarzschild black holes 
with the Gauss-Bonnet correction in the large damping limit and show that standard 
analytic techniques cannot be applied in a straightforward manner to the case of 
infinite damping. However, by using a combination of analytic and numeric tech- 
niques we are able to calculate the quasinormal mode frequencies in a range where 
the damping is large but finite. We show that for this damping region the famous 
ln(3) appears in the real part of the quasinormal mode frequency. In our calculations, 
the Gauss-Bonnet coupling, a, is taken to be much smaller than the parameter 
which is related to the black hole mass. 



1. Introduction 



Although the highly damped quasinormal modes (QNMs) of black holes are not directly 
observable, these modes have been the subject of intense investigation since the observation 
[H [2] that they may in principle provide information about the semi-classical quantum 
spectrum of black holes. This observation was motivated in large part by the special value 
of the QNM frequencies of Schwarzschild black holes in the high damping regime, namely: 

un^-^ H3)nf,/h, (1.1) 

where uj is the damping and Thh is the black hole temperature. We refer the reader to the 
literature for details [3]. The all important ln(3) in the coefficient is present for a large 
class of single horizon black holes [H Ej El [7j but not for multi-horizon black holes. Since 
the arguments that relate the QNM frequencies to black hole spectra are strictly speaking 
only valid as given in the single horizon case, it is of great interest to ask whether the 
ln(3) persists for all possible single horizon black holes of physical interest. This provides 
the motivation for the following analysis of Schwarzschild black holes with Gauss-Bonnet 
(G-B) corrections. In the rest of this paper, we will refer to these black holes as G-B black 
holes. 

G-B black holes [SI El [ID] are black holes in spacetime dimensions greater than four 
which appear when one includes, in the Einstein-Hilbert action, the leading order higher 
curvature terms arising in the low energy limit of string theories |llj. These black holes 
have recently attracted considerable attention due to the possibility of their production at 
the Large Hadron Collider if the Planck scale is of the order of one TeV [12] as occurs in 
extra dimensional brane-world scenarios. 

Reference [T3j for the first time discusses the scalar QNMs in the background of G-B 
black holes and evaluates these modes for a few values of the G-B coupling in five and 
six dimensions. The scalar QNMs in the background of G-B and charged G-B black holes 
have been studied in [T3| and this work has been extended to include G-B-de Sitter and 
G-B-Anti-de Sitter backgrounds in \T5\. The QNMs for tensor perturbations of G-B black 
holes in 5, 7, and 8 dimensions were calculated in [16J using a third order WKB formalism. 

More recently, there has been an attempt by Chakrabarti and Gupta [17] to obtain an 
analytic expression for the highly damped QNM frequencies of G-B black holes. These 
authors calculated the highly damped frequencies of tensorial perturbations in the limit 
where ^ <C 1 in which a is the G-B coupling constant and r is the radial coordinate. 
This approximation, however, is not valid in the highly damped limit where we need to 
extrapolate the WKB solutions to a region where r ~ 0. The fiaw in the results of [T7] has 
also been pointed out by Moura and Schiappa [T8] . 

The purpose of the present paper is to analyze the QNMs of D-dimensional G-B black 
holes in the large damping limit. We will show that standard analytic techniques cannot 
be applied in a straightforward manner to the case of infinite damping. However, by using 
a combination of analytic and numeric techniques we calculate the QNM frequencies in a 
range where the damping is large but finite. 

In Sec. 2, we introduce briefiy the G-B black hole and describe the general formalism 
for calculating QNMs. The specific equations describing the vector perturbations of G-B 
black holes are presented in Sec. 3, where we argue that the standard analytic techniques 
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cannot extract the infinitely damped modes. In Sec. 4 we calculate the corresponding 
QNMs in the "intermediate" damping region. The last section closes with conclusions and 
prospects for further work. 



2. Basic Formalism 

The Einstein-Hilbert action in the presence of the G-B term has the form 



a 



(D-3){D-4) 



(2.1) 

where D > 5 and a is the G-B coupling constant which is positive valued and is related 
to the Regge slope parameter or string scale. The metric for the spherically symmetric, 
asymptotically flat black hole solution of the action (12.11) is [19] 



ds^ = -f[r)dt' + $i^ + r^dn'' 



fir) 



D~2 5 



where 



/(r) = 1 + Wl + 

■'^ ' 2a 2a V 



The ADM mass, M, of the black hole is related to the parameter /i by 



M 



where is the area of a unit n-sphere. 



(2.2) 



(2.3) 



(2.4) 
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The QNM perturbations of black holes are governed generically by a Schrodinger wave- 
like equation of the form 

' V{r)\ ^ = , (2.6) 



dz^ 



+ 



for perturbations that depend on time as e"*"^*. The potential V depends in general on 
the black hole solution as well as the nature of the perturbation. It is useful to define the 
tortoise coordinate z by 

dv 

dz=^ (2.7) 

For the cases of interest here, the effective potential vanishes at both the event horizon 
[z — >■ —oo) and spatial infinity [z — >■ cxd), so that the asymptotic behavior of the QNM 
solutions is chosen to be 



-oo (r Th) 



as z — > oo r — > oo 



:2.8) 
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which represents an outgoing wave at infinity and an ingoing wave at the horizon. Since the 
tortoise coordinate is multi-valued, it is more convenient to work in the complex r-plane. 
After rescaling the wavefunction = ^P/vT obtain 

^ + R{r)^ = 0, (2.9) 

where 



with 



P 2/ 4V/ 
Here prime denotes differentiation with respect to r. 



r(n 



3. Vector perturbations and the infinite damping limit 

It was shown by Dotti and Gleiser in [20], [21], and [22] that various classes of static, 
spherically symmetric black hole metric (tensor, vector, and scalar) perturbations in the 
presence of G-B term in a spacetime with dimension D > 4 are governed generically by 
a Schrodinger wave-like equation of the form (12.61) . In this paper we focus on vector per- 
turbations because, tensor and scalar perturbations cannot be done using the numerical 
techniques employed in this paper: a crucial term in the effective potential vanishes pre- 
cisely in the intermediate damping region studied in the next section. 

The effective potential for vector perturbations is given by 



Vv{r) = q^ + 



,dK.u 
dr 



(3.1) 



with Ky{r) and gt,(r) being given by 



KJr) 



2ar2 d 
D-3dr 
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(3.2) 



and 



_ (fcg-D + 3)/ [ 2a^gef (l - f? + 2«§Efr^(l - /) + 

^^y) ~ ^2 I 



[2a(l-/)+r2]2 

where A;^ = /(/ n - 1) - 1 with / = 2, 3, 4, ■ ■ ■ . 

We now consider the QNMs in the infinite damping limit where 



(3.3) 



[Im uj] 



oo . 



(3.4) 



Since we are using complex analytic techniques, in principle the behavior of U{r) on the 
entire complex plane may be relevant. However, in the infinite damping limit case, the 
u"^/ p term in R{r) will dominate U{r) everywhere, unless one of the terms in U (r) diverges. 
This can only happen at the origin, the event horizon, the extra fictitious (i.e. complex) 
horizons, and in the case of G-B black holes at the poles where r^~^ + Sa/x = and 
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[D — 3)r'^^^ + Aafi{D — 5) = 0. Note that the poles at the points where {D — 3)r^~^ + 
4a/x(D — 5) = only exist in spacetime dimensions greater than five. Since uj'^/P also 
diverges at the event horizon and extra fictitious horizons, it will dominate there as well 
in the large damping limit. Thus, only the dominant terms of U{r) near the origin and 
the poles where r^~^ + Sa/i = and {D — 3)r^~^ + 4:a^{D — 5) = are relevant in this 
limit. Therefore, for example in five spacetime dimensions, R{r) can be approximated on 
the entire complex plane by 

, u"^ 35 48a/i 

In higher spacetime dimensions this expression becomes more complicated. 
In the WKB analysis, the two solutions to Eq. (12.91) are approximated by 

= Q-i/2(x) exp [+i Q{x')dx'] , 

(3.6) 

^W(x) = Q-'/\x) exp [-i Q{x')dx'] , 
where 

Q=(,) = (^r^^ (3.7) 

is shifted by l/(4r^) in order to guarantee the correct behavior of the WKB solutions at 
the origin. In Eq. (13. 6p . t is a simple zero of Q^. 

In order to find the QNM frequency, u, first we need to determine the zeros and poles 
of the function Q and the subsequent behavior of the Stokes and anti-Stokes lines in the 
complex r-plane. Stokes lines are the lines on which the WKB phase (/ Qdr) is purely 
imaginary and anti-Stokes lines are the lines on which the WKB phase is purely real. Due 
to the appearance of extra poles in G-B corrected Schwarzschild spacetimes, the topology of 
the Stokes/anti- Stokes lines in the complex plane will look very different from the topology 
of these lines in the Schwrzschild problem. A schematic behavior of Stokes/anti-Stokes lines 
for the G-B problem in the infinite damping limit is shown in Fig. [T] for five spacetime 
dimensions. Although this figure was drawn for a small compared to the horizon scale, the 
topology does not change as a is increased. The extra poles merely move outward in the 
complex plane as illustrated in Fig. [1] 

Once we determine the structure of Stokes/anti-Stokes lines, we can follow the ana- 
lytic method used in [23] to extract the WKB condition which consequently determines 
the QNM frequencies in the highly damped limit. In this method one needs to find a 
contour that starts and ends on an anti-Stokes line that extends to infinity in the complex 
plane, encircling only the horizon. The solution at infinity is fixed due to the boundary 
condition at infinity and the monodromy of the solution along the entire contour can be 
calculated using standard ("Stokes phenomena") rules, provided the contour stays along 
anti-Stokes lines everywhere except perhaps near zeros of Q^. This monodromy must equal 
the monodromy close to the horizon, which is in effect fixed by the boundary condition 
near the horizon. From this equality one gets an algebraic condition that determines the 
QNM frequency spectrum. One potential loop of this type is shown in Fig. [H 

Unfortunately, it seems impossible to extract a WKB condition on the highly damped 
QNM frequency given the topology shown in Fig. [H The reason is that any closed contour 
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Figure 1: A schematic illustration of the Stokes (dashed) and anti-Stokes (solid) lines in the 
complex r-plane for the G-B problem in five spacetime dimensions in the infinite damping 
limit. The open circles represent the zeros, while the filled circles are the poles of . 

beginning and ending at infinity, and circling only the horizon necessarily passes through 
the extra poles which appear due to the G-B correction. This happens for example when 
we move along the anti-Stokes lines from to t/^ and also from to ^2- At these poles the 
amplitude in the WKB solutions (13.61) goes to zero which means that after going through 
these poles we have to start with an arbitrary linear combination of the WKB solutions 
(13. 6p . In other words we lose the information imposed by the boundary condition at infinity 
as soon as we cross the poles. We are not aware of any other method which can make it 
possible to solve for the highly damped QNM frequency in the Schwarzschild spacetime 
with G-B correction. 

4. Vector perturbations in the intermediate damping region 

In this section, we investigate the approximation suggested in [17], namely 

/(r)-l-^ + ^. (4.1) 

This approximation is only valid if all WKB calculations are done in a region of the complex 
plane where ^ 1 or equivalently r ^ (8a/i)^/*^^~^\ Thus, all the relevant Stokes/anti- 
Stokes lines and the zeros of the function Q along the contour in the WKB analysis should 
be far from the pole at the origin and at r^~^ + 8a/i = and (D — 3)r^~^-|-4a/i(L) — 5) = 0. 
This implies that the approximation (14. ip . forbids taking the infinitely damped limit of 
the QNM frequency because in this limit the zeros of the function Q approach the origin 
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Figure 2: The zeros of the function Q in five spacetime dimensions for different values of 
\uj\. 



where the approximation (14. ip becomes invahd. The authors of [T7j ignored this crucial 
point in their calculations and consequently arrived at incorrect conclusions. 

However, if we lower the damping, the zeros of the function Q will move away from the 
poles. Fig. [2] shows how the zeros of the function Q in five spacetime dimensions move in 
the complex plane as the damping decreases. In this figure 2/i = 1, a = 10~^, and 1 = 2. 
As one can see, six of the zeros of the function Q will move significantly outward away 
from the poles in the complex plane. Four of the zeros will stay relatively close to the poles 
at + 8a/i = 0. Thus, for intermediate damping, the six zeros move to a region where 
r ^ (8a/i)^/^ and the approximation fl4.ip is valid. It is easy to show that, in general, for 
the approximation (14. ip to be valid we need to have 



\UJ 



[Im cuf < 



(4n2 - 
(8a/i)"+i 



(4.2) 



Note that as in the purely Schwarzschild calculation, we also require that the damping to 
be large compared to the real part of the frequency, which implies that 



UJ » 



(4.3) 



and guarantees that the zeros along the contour used are well inside the horizon radius. 

In this intermediate damping region the topology of the Stokes/anti- Stokes lines will 
look very different from the one in Fig. [TJ Fig. [3] shows the numerically generated topology 
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of the Stokes/anti-Stokes lines for five dimensional G-B black holes for an intermediate 
damping of |ci;| = 10^. In this figure 2/i = 1, a = 10~^, and 1 = 2. To generate this figure 
we have written a numerical recipe that uses the approximation 



a+dz 



2 1 



(4.4) 



when the stepsize \dz\ is small enough. Here Q = Q^. The numerical calculation starts at 
the zeros of the function Q. To find the Stokes lines, the program finds the point where 



iQia)+Q'{a)dzY/^-Qiaf' 



3Q'{a) 

and to find the anti-Stokes lines the program finds the points where 



< e 



(4.5) 



T [2 1 
e < Im < - — 



3g'(a) 



{Q{a) + Q'{a)dzf'-Q{af/^ 



< e . 



(4.6) 



Here e ideally needs to be zero, but in the numerical calculations we take it to be a small 
number which generally works well if we take it to be roughly equal to the stepsize \dz\. 
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Figure 4: A schematic illustration of the Stokes (dashed) and anti-Stokes (solid) lines in 
the complex r-plane for the G-B problem in five spacetime dimensions for an intermediate 
damping. The open circles represent the zeros, while the filled circles are the poles of . 

We plot the schematic behavior of the Stokes/anti-Stokes lines in Fig. IHfor G-B black 
holes in five spacetime dimensions in the intermediate damping region. Note that the 
topology of the Stokes/anti-Stokes lines in this damping region looks almost identical to 
the topology of these lines for five dimensional Schwarzschild black holes. The only region of 
the complex plane which is different from the Schwarzschild problem is the region contained 
inside the anti-Stokes lines connecting the outer six zeros. This region stays outside of the 
loop, shown in Fig. HJ that we follow to extract the WKB condition and therefore will 
not affect the QNM calculations. In addition, compared to the Schwarzschild problem, 
there exist four extra Stokes lines in Fig. H] that extend to infinity. These extra Stokes 
lines, which cross the anti-Stokes lines, will not affect our calculations either, because the 
exponential dominance of the WKB solutions (13. 6p only changes at the so called "transition 
points" which are the zeros of the function Q for this particular problem. More specifically, 
along the anti-Stokes line connecting t\ to ^3 there will be no change in the dominance of 
the WKB solutions due to the crossing of the Stokes line. 

Let us see how the function R is modified if we use the approximate metric function 
/(r) in Eq. (14.11) . In the large but finite damping limit and when <C 1, the function 
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Figure 5: 2??-^^^= Re a; as a function of — Im uj for vector perturbations in five spacet 
dimensions. 
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R takes the form 



R{r) = J, U 



{An" - l)/i^ 



1 + 



n 



9n^ + 18n + 10 fSafi 



2(4^2 - l)(n - 1) Vr"+ 



+ 



(4.7) 

where n = D — 2. It is easy to see that as long as the calculations are done in a region of 
the complex plane where r ^ (Sa/i)^/"^^, we get the same function R as in Schwarzschild 
black holes plus the leading terms which add a small correction due to the presence of the 
G-B coupling constant a. 

In Fig. we plot the QNM frequency as a function of damping. Note that according to 
this figure the deviation of the real part of the frequency from ln(3) due to G-B corrections 
is very small. To generate the numerical data in this figure, we have followed the same 
combination of analytic and numeric techniques developed in [23]. First, we use the analytic 
method of ^23j to find the general WKB condition 



_g2i743 _ 



;i + e2*T«) 

p2i7i3 



for the Schwarzschild-type topology shown in Fig. HI Here 



7a6 



ta 



Q{r)dr 



(4.9) 



and F is the integral of the function Q along a contour encircling the pole at the event 
horizon in the negative direction. Using the residue theorem we can evaluate F in five 
spacetime dimensions: 



Qdy 



-2TTi Res Q 

r=\rh\ 



-2m 



2fi + a 
2V2/i - a 



(4.10) 
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The phase integrals 713 and 743 are calculated numerically using the approximation (14.41) . 
For the particular case in Fig. [3] we find 713 ~ —3.12853 ■ ■ ■ radians and 743 ~ —3.16878 ■ • • 
radians. Therefore we get 

g2*r ^ ^^^^T^^ ^ -2.99494 ■ ■ ■ + iO. 16088 ■ • • . (4.11) 

In other words 

27T ^^^^^ UJ ^ ln(2.99926 ■ • •) - i{2n + l)7r , (4.12) 
-\/2yU — a 

where n ^ 1. In Fig. [3l we have taken \uj\ ~ — Im u = 10^ which means that n ^ 10^. 



5. Conclusion 

We have calculated, using a combination of analytic and numeric techniques, the QNM 
frequencies for G-B black holes in the intermediate damping region, in which the imaginary 
part of the QNM frequency is large relative to the black hole mass, but not infinitely large. 
Physically this corresponds to taking the low energy limit, a ^ before u; ^ 00. As 
expected in this limit, the frequencies are perturbations of the ln(3) that appears for the 
highly damped QNM perturbations of the Schwarzschild black hole. Interestingly, the G-B 
corrections decrease the real part of the frequency as the damping is increased. Therefore, 
contrary to what was claimed by the authors of [17] , we conclude that the "universal" value 
of ln(3) with a small deviation due to G-B corrections does appear in G-B black holes. 

The story will change for the infinite damping limit. The infinitely damped QNMs 
essentially "see" the full structure of the solution including the extra poles of the function Q 
which only exist in G-B black holes. These extra poles will affect the value of the frequency 
in this limit. This phenomenon is very similar to what happens in Reissner- Nordstrom (R- 
N) black holes with small charge, where the real part of the QNM frequencies makes a 
transition from R-N value in the infinitely damped limit to the Schwarzschild value of 
ln{3)Thh for an intermediate damping [2lj. In that case also the infinitely damped QNMs 
essentially "see" the full structure of the solution, including the inner horizon, which affects 
the value of the frequency in this limit. The intermediate damping on the other hand only 
probes the outer horizon in R-N black hole, and for this range the "universal" value of 
ln(3) is reproduced. In the present context, this behavior suggests that in order to solve 
for the highly damped QNMs we may need to work with a more complete model of black 
holes in string theories rather than the G-B black hole which only includes the leading 
order higher curvature terms arising in the low energy limit of string theories. 
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